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Abstract 

Both scalar and Dirac quasinormal modes in Garfinkle-Horowitz-Strominger black hole space- 
time are studied by using the WKB approximation and the Poschl- Teller approximation. For scalar 
field, we find that the QNMs with higher dilatons decay more rapidly than that with lower ones. 
However, this is not the case for mass m. Fields with higher mass will decay more slowly. The 
similar behaviors appear in the case of Dirac field. We also find that QNM frequencies evaluated 
by using WKB approximation and Poschl- Teller approximation have a good agreement with each 
other when mode number n is small. 

PACS: number(s): 04.70.Dy, 04.70.Bw, 97.60.Lf 

I INTRODUCTION 

The stability of a black hole has been discussed for many years. The pioneering work on this 
problem was carried by Regge and Wheeler PE], who studied the linear perturbation of Schwarzschild 
black hole. VishveshwarajBj and Chandrasekhar^ furthered this work by bringing forward quasinormal 
modcs(QNMs). Being its importance both in the analysis of the stability of the black holes and in the 
search for black holes and its gravitational radiation, the QNMs of black holes in the framework of 
general relativity have been studied widely. 

QNMs of black holes are defined as proper solutions of the perturbation equations belonging to 
certain complex characteristic frequencies which satisfy the boundary conditions appropriate for purely 
ingoing waves at the event horizon and purely outgoing waves at infinity Generally speaking, the 
evolution of field perturbation on a black hole consists roughly of three stages[B]. The first one is an 
initial wave dependent on the initial form of the original field perturbation. The second one involves 
the damped oscillations called QNMs, which frequencies and damping times are entirely fixed by the 
structure of the background spacetime and are independent of the initial perturbation. The last stage 
is a power-law tail behavior of the waves at very late time which is caused by backscattering of the 
gravitational field. 

QNMs were firstly used to study the stability of a black hole. Then people used it to study the 
properties of black holes for the definite relations between the parameters of the black hole and its 
QNMs. Latest studies in asymptotically flat space have acquired a further attention because of the 
possible relation between the classical vibrations of black holes and various quantum aspects, which was 
proposed by relating the real part of the quasinormal frequencies to the Barbero-Immirzi parameter, 
a factor introduced by hand in order that loop quantum gravity reproduces correctly entropy of the 
black hole[3 El Ej- Moreover, QNMs also relate to superstring theorv |T7)l 111) . This is known as the 
AdS/CFT correspondence, which argues that string theory in anti-de Sitter(AdS) space is equivalent 
to conformal field theory(CFT)in one less dimension. 

* E-mail: ygshen@center.shao.ac.cn 
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We know that the Schwarzschild spacetime can describe an uncharged black hole in string theory 
very well (except in the region near the horizon) when the mass of the black hole is larger than the 
Plank mass. However, this is not the case for a charged black hole in string theory. Gibbons and 
Maeda^l first obtained its solution in this case, latter Garfinkle, Horowitz and Strominger 13 find 
an another kind of solution, which is known as Garfinkle-Horowitz-Strominger solution. It has been 
showed that there exist many essential differences between the Garfinkle-Horowitz-Strominger black 
hole and the Reissner-Nordstrom black hole. This makes it necessary for us to study the QNMs of 
Garfinkle-Horowitz-Strominger dilaton black hole, although the similar work of Reissner-Nordstrom de 
Sitter black hole has been done by Ref.|14|. The main purpose of this paper is to study the QNMs 
of Garfinkle-Horowitz-Strominger dilaton black holes by evaluating the quasinormal mode frequencies 
both of scalar and Dirac fields. 

Many methods are available for calculating the quasinormal mode frequencies. Most of them are 
numerical in nature. However, there are two methods often used, i.e., WKB potential approximation, 
which was devised and developed in |15lHB] . and Poschl- Teller (PT) potential approximation, which was 
first proposed by Poschl and Teller 17 . Both methods have their own merits. WKB approximation is 
accurate for the low- lying modes [T%|. while the Poschl- Teller approximation is simpler than the first one 
and accurate enough for our purpose. Therefore, we utilize the first method to calculate the frequencies 
of the scalar field, but use the second method to calculate the frequencies of the massive Dirac field. 
Moreover, we use these two methods to calculate the massless Dirac field so that we can check the 
consistency between these two kinds of methods. 

The paper is organized as follows: In the next section we study the scalar QNMs in Garfinkle- 
Horowitz-Strominger dilaton black hole, and evaluate their quasinormal mode frequencies. In Sec. Ill, 
we discuss the Dirac QNMs and calculate their characteristic frequencies. Some conclusions are drawn 
in the last section. 

II SCALAR QUASINORMAL MODES IN THE 
GARFINKLE-HOROWITZ-STROMINGER BLACK HOLE 

We consider the QNMs in the Garfinkle-Horowitz-Strominger black hole described by the metric [T3*l 

ds 2 = -e 2U dt 2 + e- 2U dr 2 + R 2 (dO 2 + sin 2 9dtp 2 ) , (1) 

e - 2 = e -a*o (1 _£) (2 ) 
r 

F = Q sin Odd A dtp, (3) 

with 

= 1-^, (4) 
r 



R = y/r(r-a), (5) 

where M is the mass of the black hole, Q is magnetic charge, cj> and F are scalar and magnetic field, 
respectively. The parameter a is defined by 

a = -^-e" 200 (6) 

and 4>o is an arbitrary constant. 

The best way to deal with the QNMs for the scalar field is to solve equations deduced from 
the action. However, it is very difficult, if possible, to solve these equations because we must find 
general solutions rather than static, spherically symmetric solutions which obtained in Rcf. 13 when 
we discuss scalar QNMs in Garfinkle-Horowitz-Strominger black hole. Under this consideration, we use 
another effective way to solve this question, the way has been used to discuss the entropy in Garfinkle- 
Horowitz-Strominger black hole in Ref.|19|. 
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The scalar field with mass to satisfies the wave equation 

(□ - m 2 ) $ = 0. 

If we introduce a tortoise coordinate 



dr* 



-2U 



dr, 



and define 



Then the wave equation can be written as 



*J + (^-v)*-a 



where 



V 



2M 
r 



M(2r 



r 3 (r — a) 



2M 
r 



a 2 1(1 + 1) 

4r 2 (r - a) 2 r (r - a) 



(7) 

(8) 
(9) 

(10) 
(11) 



It is an extremal black hole for Garfinkle-Horowitz-Strominger black hole when a = 2M, since 
there exists a curvature singularity at r = a as showed in Eqs.(T5). In this case, many new unknown 
properties occur and thus further work is needed. In this article, we only calculate the quasinormal 
mode frequencies when a ^ 2M. It is obvious that the effective potential satisfies 



V 
V 



e 2M . 

™2 



as 
as 



-oo, 

-CO. 



As to is increased, the peak value of the effective potential will be smaller than the asymptotic value 
to 2 . Reference |20|pointed out that there will be no quasinormal modes if the peak of the potential 
is lower than the asymptotic value m 2 , i.e., there exists a maximum value above which quasinormal 
modes can not occur. Therefore, one can estimate its value m max by using the relation 



V(r 



max i " L max 



,l,a) 



(12) 



The results are tabulated in Table I. Here we have used the mass M of the black hole as a unit of mass. 
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Table I: Maximum values of the mass m of 
the scalar field above which QNMs can not 
occur. 



where 



From the discussion above 
must let to smaller than m max 
the quasinormal mode frequencies 



we see that we 
when we calculate 
One can obtain 



the complex quasinormal mode frequencies uj by us- 
ing WKB approximation. The formula, carried to 
third order beyond the eikonal approximation, is given 
by 16 



A = 



y + (-2F ") 3 A 




\n-\)(-2Vi;\- O). 



(7 + 60a z 



(13) 



(14) 
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Table II: Quasinormal mode frequencies for scalar field 
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Figure 1: Scalar quasinormal mode frequen- 
cies 
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(n) 
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(16) 
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r,=r»(r ) 

We can obtain the complex quasinormal mode frequencies by plugging the potential V in eq. Ijll|l 
into the formula above. The values for < I are listed in Table II. The values for negative n are related 
to those with positive n by reflection of the imaginary axis. Figure 1 is its corresponding figure. 



Ill DIRAC QUASINORMAL MODES IN THE 
GARFINKLE-HOROWITZ-STROMINGER BLACK HOLE 

The Dirac equation in a general background spacetime is|2] 

[ 7 a < (d M + r„) + m] *■ = o, 

where is the inverse of the tetrad defined by the metric g^ v , 

and is the spin connection which is defined as T M = g [7 a ,7 fc ] e^a&bv.^ ■ 
If we define 



-u/2 ( ^G ± (r)^(e, ^ p)/R 



and use the tortoise coordinate dr* 



-2U 



F ± (r)^(e, V )/R )" 
dr, the Dirac equation can be decoupled as 



d 

dr* 
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G 



k/R 



-k/R 



where k satisfies 



F 
G 



3 = I + \ 

3 = l-h 
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UJ 
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If we make a change of variables 



where = tan" 1 (^) 
Eq.(21) becomes 



F 
G 



d 
dr* 



(!) 



cos 
— sin 



(f) 



F 
G 



F 
G 



W 



-F 
G 



G 

—F 



(18) 
(19) 

(20) 

(21) 

(22) 
(23) 

(24) 
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where we have made a coordinate transformation 

. 1 

2uj 



-tan 



mR 
K 



and 



W 



( K /Rf 



1 + nmR'e 2U / [2uj (k 2 + m 2 R 2 )} ' 



Eq.(24) can also be written as 



d 2 F 

df 2 

d 2 G 
df 2 



f (lo 2 -Vi)F = 0, 
+ (lu 2 -V 2 )G = 0, 



and 



Vi, 2 = ± 



dW 



(25) 
(26) 

(27) 
(28) 

(29) 



where V1.2 are supersymmetric partners derived from the same superpotential. It is obvious that 
potentials related in this way possess the same characteristic frequencies for both F and Gj^i which 
means that Dirac particles and antiparticles have the same quasinormal mode spectra in the Garfinkle- 
Horowitz-Strominger black hole spacetime. Therefore, it is reasonable to concentrate just on Eq.(27) 
in evaluating the quasinormal mode frequencies in the next sections. 



A Quasinormal mode frequencies for the massless Dirac field 



In this subsection, we shall evaluate the quasinormal mode frequencies of the massless Dirac field 
in the Garfinkle-Horowitz-Strominger black hole as a ^ 2M by using WKB potential approxima- 
tion and Poschl- Teller potential approximation, so that we can compare the results with each other. 
Zhidenko 22 has proved that the quasinormal mode frequencies obtained by the Poschl- Teller approx- 
imation and the sixth-order WKB approximation are in a very good agreement for the Schwarzschild 
de Sitter black hole. 



We rewrite Eq.(27) as 



where 



d 2 F 
drl 



V 



R 2 



V)F = 0, 



±[Re u U ~e u R 



(30) 



(31) 



Here we have written V± as V because we shall not consider V% , which will give the same spectrum 
of quasinormal mode frequencies for the previous reason. 

If we let a = 0, the equations 1)30(1 - Q31JI give the results of the Schwarzschild black hole 20 . In the 
following steps, we will use the equations l(30(l - ((31|l to evaluate the quasinormal mode frequencies when 
a ^ 2M. 

From Fig. 2, we know that the effective potential V, which depends only on the value of r for fixed 
K I and a, has a maximum over r G (2M, +oo). The location ro of the maximum has to be found 
numerically. But for | k |— ► 00, we can get the position of the potential peak from V1.2, 



ro (I « 



00 



a + 6M 




and 



(I K I) < r (I K H °°) 

ro (| K |) > r (| k |-> 00) 



32aAf 
(a + 6Mf 

(for Vi), 
(for V 2 ). 



(32) 



G 



If we treat the mass M of the black hole as a unit of mass and length, the value of the ro (| k — > oo) 
only depends on the coupling coefficient a. It gives ro (| k |— ► oo) = 3 for a black hole with a = 0, 
which agrees with the result for the Schwarzschild black hole 20 . From Fig. 3, we see that the peak 
value of the effective potential V increases with a. 
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Figure 2: Variation of the effective poten- 
tial V for massless Dirac field V with n as 
a= 1. 



Figure 3: Variation of the effective poten- 
tial V for massless Dirac field with a as 

K= 1. 



1 Dirac quasinormal mode frequencies by using WKB approximation 

We now evaluate the quasinormal mode frequencies for the massless Dirac field in Garfinkle- 
Horowitz-Strominger black hole using the WKB potential approximation. In Sec. II, we have used this 
kind of approximation to evaluate the quasinormal mode frequencies for scalar field. We just need plug 
the effective potential V(r, a, k) in Eq.(31) into Eq.(17) and use Eqs.(13)-(16) to evaluate the complex 
frequencies for Dirac field. Here we only list the results for positive n in Table V, the corresponding 
values for positive k are plotted in Fig. 8. 



2 Dirac quasinormal mode frequencies by using Poschl- Teller approximation 

The potential calculated by Eq.(29) are smooth function, integrable over the range of r* G 
(— oo, +oo), i.e., 



V 
V 



e 2M . 

_-2 



as 

as r* 



-oo, 
-oo. 



From this relation, we know that the effective potential drops exponentially to zero as r* 



(but falls off as r as 



Poschl- Teller approximation potential, |17l l2"3*| 

Vpt — 



+oo). So we can study the Dirac quasinormal modes by using the 



cosh a(r* — ?*o) 



(33) 
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The quantities Vq and a > are given by the height and curvature of the potential at its maximum 
(r* = r ), i.e., 



V = V(r ), 



a 



1 



dri 



(34) 



The quasinormal mode frequencies of the Poschl- Teller potential can be evaluated analytically 23 



w„ = \j V - — - a (n 



1/2)*, 



(n = 0,l,2,--.). 



(35) 



The results for positive n are listed in Table VI. From Table V and Table VI, we see that the 
values evaluated by this two kinds of approximation have a good agreement with each other. The real 
part of the complex frequencies, which decreases slowly as mode number n increases, keeps unchanged 
as n changes as showed in Fig. 8. However, as to the magnitude of the imaginary part of the frequencies 
evaluated by these two kinds of approximations, both increase with n. Therefore, both approximations 
will break down when n is large. In this article, we see that the accuracy of two methods is good 
enough for any n < k. It is also showed in Fig. 8 that the complex frequencies increases with a. This 
indicates that the QNMs with higher dilaton decay more rapidly. 



B Quasinormal mode frequencies for the massive Dirac field 

In this subsection, we shall evaluate the quasinormal mode frequencies of the massive Dirac 
field as a / 2M in the Garfinkle-Horowitz-Strominger black hole by using Poschl- Teller potential 
approximation. 



We can rewrite Eq.(27) as 



d 2 F 
dri 



where 



(lo 2 -V)F = 0, (36) 
2 In (r/2 - 1) + i tan" 1 (my/r(r - o)/k) , (37) 



p 2U ( K 2 , 2p2y/* 

V = 6 {K +m R) f{e u U'R( K 2 + m 2 R 2 ) 

R 2 [(K 2 + m 2 R 2 ) + e 2U mKR , /2cj} 2 1 V ' 

- K 2 P U W - 1 e 3U nmR r, 2 2 R 2\ 

2uj{k 2 + m 2 R 2 ) + e 2 U Km R'/2u; [[ + ) 

■ (2U'R' + R") - 2m 2 RR' 2 ] + (k 2 + m 2 R 2 f /2 } . (38) 

The dependence of potential V(r, to, k, a, u>) on m, k, a and ui is showed in Fig. (4-7), respectively. 
From Fig. 4, we can see that the peak of the potential increases with to. However, the potential has an 
asymptotic value m 2 as r — > 00. Hence, there exists a maximum m max above which the height of the 
peak is lower than the asymptotic value to 2 . As a result, there will be no quasinormal modes when to 
is larger than m max . So we can use the same method which we have used in Sec. II to estimate this 
maximum value. 

V (^rnaxj TTlmaxi K^CL^U] — ^mor) — y^max) ■ (*^) 

The values calculated in this way are tabulated in Table III and Table IV, where we have also 
given the maximum values of fi — m/n which will be used in the following discussions. 

Fig. 5 shows the dependence of the potential on the angular momentum quantum number k, from 
which we see that the behaviors of the potential in this case are similar to that of the massless one as 
showed in Fig. 2. The same case happens in the relation between potential V and a, as showed in Fig. 3 
and Fig. 6. The variation of the potential with a for massive Dirac field approaches to that of massless 
one. However, Fig. 7 shows that the potential for massive Dirac field depends obviously on the energy 
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Figure 4: Variation of the effective poten- 
tial V for massive Dirac field with m as 
k = l,a = l,u> = 1. 



Figure 5: Variation of the effective poten- 
tial V for massive Dirac field with k as 
m = 0.2, a = l,u> = 1. 



(or frequency) to. This tells us that we can't evaluate the quasinormal mode frequencies as done for 
the massless ones. In the following parts, we'll evaluate the complex quasinormal mode frequencies for 
massive Dirac field by using Poschl- Teller potential approximation. 

From Table III and IV, we can see that | /i maa; |< 1 for any k ^ 0, and as | n |— > oo, /x — » 
regardless of the value of m. Therefore, we can evaluate the quasinormal mode frequencies in the 
Poschl- Teller approximation as power series of /i for given values of K as showed in [201 1241 12*5] . The 
position r max of the peak of the massive potential can be obtained by expanding r as power series of 
/i up to order /i e , 

Trnax = r a + r X fJL + r 2 ^ 2 + r 3 ^ 3 + r 4( u 4 + r 5 /i 5 + r 6 ^ 6 = r + E, (40) 

and requires 

= V (w) = V (r ) + EV" (r ) + isV" (r ) + ^ 3 V W (r ) + 

2l E4y(5) (ro) + ^ 5V{6) (ro) + W^ y(7) (ro) ' (41) 
where ro is the position of the peak value for the massless case obtained in the last section. The 
coefficients can be then evaluated order by order from this equation. In the following step, we 
expand lo up to order /i 6 

10 = LUq + LOXfJ* + L0 2 ^ 2 + L03fJ? + LO^fl + L0 5 [1 5 + L0 6 [1 6 , (42) 

and plug it into r max , then expand it again up to order ^t 6 . Using this new expansion for r max , we can 
expand the right hand side of Eq.(35). Note that Eq.(34) should be written as 



V = V{r max ), a 2 = --i r 

ZVq 



Then we can solve u>i self-consistently order by order in fi. The values evaluated in this way are listed in 
Tables (VII- VIII). Quasinormal mode frequencies for any fixed /i < fi ma x can be obtained by plugging 
these coefficients into Eq.(42). Fig. 9 is the corresponding figure, from which, we can see that the real 
parts of the quasinormal mode frequencies increase with m, but the magnitude of the imaginary parts 
of the frequencies decrease as m increase. This means that field with higher mass will decay more 
slowly. 
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4 


0.9436 


0.2359 


1.0158 


0.2539 


1.1165 


0.2791 


5 


1.1922 


0.2384 


1.2823 


0.2565 


1.4080 


0.2816 



Table III: Maximum values of the mass m and /x(= m/n) of the Dirac field above which QNMs can't 
occur for positive k. 





a 


= 


a = 


0.5 


a 


= 1 


K 


m 




m 




m 




1 


0.3333 


-0.3333 


0.3524 


-0.3524 


0.3780 


-0.3780 


2 


0.5723 


-0.2862 


0.6089 


-0.3045 


0.6589 


-0.3295 


3 


0.8190 


-0.2730 


0.8732 


-0.2911 


0.9479 


-0.3160 


4 


1.0673 


-0.2668 


1.1394 


-0.2849 


1.2389 


-0.3097 


5 


1.3164 


-0.2633 


1.4063 


-0.2813 


1.5306 


-0.3061 



Table IV: Maximum values of the mass m and /i(= m/n) of the Dirac field above which QNMs can't 
occur for negative n. 
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IV SUMMARIES AND DISCUSSIONS 



Both scalar and Dirac quasinormal mode frequecies in Garfinkle-Horowitz-Strominger black hole 
are evaluated by using WKB approximation and Poschl- Teller approximation. As for scalar field, the 
complex frequencies depend on the mass m of the field, the orbital angular momentum I, and the 
coupling coefficient a. From the results we have obtained in Table II and Fig.l, we can see that the 
real parts of the frequencies increase with m and a, but they decrease with n. However, the magnitude 
of the imaginary parts of the frequencies increase with a and n, but decrease with m. This indicates 
that field with higher mass will decay more slowly, while QNMs with higher dilaton will decay more 
rapidly. The similar behaviors appear in the Dirac field. But for massive Dirac cases, we have made 
an expansion in powers of the parameter [i = m/n. In this way, we can obtain the numerical values of 
the quasinormal mode frequencies. 

Another conclusion is that the quasinormal mode frequencies evaluated by using WKB approxi- 
mation agree with the results obtained by using Poschl- Teller approximation when mode number n is 
small, as showed in Fig. 8. It is important that we have the same conclusion for massive cases. 
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Tabic V: Dirac quasinormal mode frequencies by using WKB approximation 



a = 



a = 0.5 



a = 1 



a = 1.5 



0.1765-0. lOOli 
0.3786-0. 0965i 
0.3536-0. 2988i 
0.9602-0.0963i 
0.9496-0. 2902i 
0.9300-0.4876i 
0.9036-0. 6892i 
0.8721-0. 8944i 



0.1976 

0.4168- 

0.3947- 

1.0539- 

1.0445- 

1.0269 

1.0032 

0.9749- 



0.10251 
0.0994i 
0.3061i 
0.0992i 
0.2987i 
0.5015i 
0.7081i 
0.91801 



0.2275-0. 1054i 
0.4716-0. 1028i 
0.4534-0.3147i 
1. 1889-0. 1025i 
1.1810-0.3088i 
1. 1661-0. 5175i 
1. 1460-0. 7297i 
1.1222-0.945H 



.2780-0. 1083i 
5651-0. 1061i 
.5536-0.3226i 
4192-0. 1058i 
.4137-0.3180i 
4037-0.5323i 
.3902-0.7493i 
.3744-0.9690i 



Table VI: Dirac quasinormal mode frequencies by using Poschl- Teller approximation 



K 


n 


a = 


a = 0.5 


a = 1 


a = 1.5 


1 





0. 1890-0. 1048i 


0.2083-0.1064i 


0.2360-0. 1082i 


0.2829-0. 1092i 


2 





0.3855-0.0991i 


0.4229-0.1015i 


0.4766-0. 1043i 


0.5684-0. 1068i 




1 


0.3855-0. 2972i 


0.4229-0.3045i 


0.4766-0.3130i 


0.5684-0.3203i 


5 





0.9625-0.0966i 


1.0561-0.0995i 


1. 1908-0. 1028i 


1.4204-0. 1059i 




1 


0.9625-0. 2899i 


1.0561-0.2984i 


1.1908-0.3084i 


1.4204-0.3176i 




2 


0.9625-0.4832i 


1.0561-0.4974i 


1.1908-0.5140i 


1.4204-0. 5293i 




3 


0.9625-0. 6764i 


1.0561-0.6964i 


1. 1908-0. 7195i 


1. 4204-0.741 li 




4 


0.9625-0. 8697i 


1.0561-0.89531 


1.1908-0.925H 


1.4204-0.9528i 




















a-Q a=0 a=0 . 5 1 15 










oaa » » » » • . » » • • 
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Figure 8:Massless Dirac quasinormal mode frequencies evaluated by using WKB and Poschl- Teller 

potential approximation 
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Tabic VII: The real parts of the quasinormal mode frequencies for massive Dirac field 



a 


K 


n 




LU 








LU 2 


L03 


a>4 


U> 5 


U> 6 





1 








1890 


-0 


1306 


1 


1736 


-1.2880 


-1.2438 


25.707 


-104.57 




2 








3855 


-0 


1583 


1 


8235 


0.3132 


-0.1220 


2.0807 


4.3001 






1 





3855 


-0 


1207 


1 


8048 


-0.4484 


2.3177 


-4.2165 


-8.5152 




5 








9625 


-0 


1654 


4 


3544 


0.4607 


2.6316 


1.8376 


14.149 






1 





9625 


-0 


1572 


4 


3535 


0.2019 


2.9479 


2.7Q2" 


6.5863 






2 





9625 


-0 


1439 


4 


3510 


-0.1557 


3.4228 


1.9923 


-0.4496 






3 





9625 


-0 


1294 


4 


3472 


-0.4621 


3.8330 


-1.0875 


-0.8916 






4 





9625 


-0 


1160 


4 


3433 


-0.6652 


4.0637 


-5.3255 


3.8645 


0.5 


1 








2083 


-0 


1350 





9330 


0.0801 


-0.2485 


-1.0873 


-0.3100 




2 








4229 


-0 


1497 


1 


5490 


0.1968 


0.2509 


1.2826 


2.5888 






1 





4229 


-0 


1166 


1 


5328 


-0.3105 


1.8138 


-2.5407 


-2.8215 




5 





1 


0561 


-0 


1555 


3 


7113 


0.2922 


2.3756 


1.1681 


9.1466 






1 


1 


0561 


-0 


1484 


3 


7104 


0.1269 


2.5702 


1.4507 


5.7704 






2 


1 


0561 


-0 


1368 


3 


7083 


-0.1086 


2.8675 


0.9400 


2.4967 






3 


1 


0561 


-0 


1236 


3 


7051 


-0.3182 


3.1369 


-0.7745 


2.1499 






4 


1 


0561 


-0 


1111 


3 


7017 


-0.4620 


3.3019 


-3.1439 


4.3852 


1 


1 








2360 


-0 


1265 





7241 


0.0224 


0.0763 


-0.4976 


0.5058 




2 








4766 


-0 


1385 


1 


2327 


0.0933 


0.4391 


0.6423 


1.5177 






1 





4766 


-0 


1110 


1 


2196 


-0.1944 


1.2855 


-1.3001 


0.1945 




5 


o 


1 


1908 


-0 


1430 


2 


9092 


0.1280 


3.0462 


1.0794 


-3.0867 






1 


1 


1908 


-0 


1439 


2 


8908 


-0.1161 


3.3237 


4.6232 


-2.5366 






2 


1 


1908 


-0 


1452 


2 


8650 


-0.5293 


3.5692 


9.5014 


-0.0203 






3 


1 


1908 


-0 


1468 


2 


8439 


-1.0123 


3.5679 


13.628 


3.5289 






4 


1 


1908 


-0 


1484 


2 


8326 


-1.4886 


3.3481 


16.245 


6.0940 



n=0 


0.2 
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0.2 ' . 
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'' 


0.4 


n=2 






0' 


0.4 


n=3 






o' 


0.4 


n=4 











0.2 0.4 0.6 0.8 1 1.2 

Re(w) 

Figure 9: Massive Dirac quasinormal mode frequencies 
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Tabic VIII: The imaginary parts of the quasinormal mode frequencies for massive Dirac field 



a 


K 


n 




LU 


LUl 




LU 2 




LU3 


LU4 


w 5 







1 





-0 


1048 


-0.0629 





4857 





7429 


-3.6364 


15.433 


41.564 




2 





-0 


0991 


-0.0292 





6419 





5097 


0.4961 


3.9247 


18.790 






1 


-0 


2972 


-0.0586 


1 


9232 





6183 


4.1854 


12.488 


3.3177 




5 





-0 


0966 


-0.0112 





7081 





2414 


1.0021 


1.0909 


9.5112 






1 


-0 


2899 


-0.0310 


2 


1245 





6171 


3.1196 


4.5254 


23.934 






2 


-0 


4832 


-0.0450 


3 


5407 





7701 


5.5208 


9.4815 


29.923 






3 


-0 


6764 


-0.0528 


4 


9555 





7419 


8.2001 


13.885 


32.585 






4 


-0 


8697 


-0.0559 


6 


3684 





6276 


11.014 


16.435 


37.898 


0.5 


1 





-0 


1064 


-0.0555 





4235 





4624 


-1.1012 


5.0715 


0.3089 




2 





-0 


1015 


-0.0262 





5267 





3303 


0.5246 


2.5273 


9.7331 






1 


-0 


3045 


-0.0547 


1 


5759 





4160 


3.2193 


7.3289 


5.7093 




5 





-0 


0995 


-0.0101 





5766 





1567 


0.9083 


0.7544 


6.2479 






1 


-0 


2984 


-0.0284 


1 


7298 





4049 


2.7927 


2.8027 


16.737 






2 


-0 


4974 


-0.0418 


2 


8828 





5128 


4.8468 


5.5232 


23.440 






3 


-0 


6964 


-0.0498 


4 


0346 





4989 


7.0743 


7.9666 


28.421 






4 


-0 


8953 


-0.0535 


5 


1849 





4213 


9.3972 


9.4492 


34.445 


1 


1 





-0 


1082 


-0.0468 





3346 





2387 


-0.4409 


2.8934 


-1.6048 




2 





-0 


1043 


-0.0225 





4041 





1797 


0.4793 


1.3415 


4.2433 






1 


-0 


3130 


-0.0495 


1 


2077 





2342 


2.2696 


3.5442 


5.2856 




5 


o 


-0 


1028 


0.0012 


o 


4498 


o 


3779 


0.7181 


-3.0411 


2.0343 






1 


-0 


3084 


0.0035 


1 


3419 


1 


0698 


2.3336 


-7.6403 


5.5251 






2 


-0 


5140 


0.0053 


2 


2172 


1 


6034 


4.2615 


-9.1513 


9.3059 






3 


-0 


7195 


0.0065 


3 


0780 


1 


9521 


6.3210 


-8.1271 


15.282 






4 


-0 


9251 


0.0071 


3 


9316 


2 


1405 


8.3013 


-5.8811 


23.284 
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